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\Q Consider the continuum of points along the edges of a network, i.e., an undirected graph with positive edge 

weights. We measure distance between these points in terms of the shortest path distance along the network, 
known as the network distance. Within this metric space, we study farthest points. 

We introduce network farthest-point diagrams, which capture how the farthest points — and the distance 

to them — change as we traverse the network. We preprocess a network G such that, when given a query point 

c/3 q on G, we can quickly determine the farthest point (s) from q in G as well as the farthest distance from q 

i in G. Furthermore, we introduce a data structure supporting queries for the parts of the network that are 

farther away from q than some threshold R > 0, where R is part of the query. 

We also introduce the minimum, eccentricity feed-link problem defined as follows. Given a network G with 
geometric edge weights and a point p that is not on G, connect p to a point q on G with a straight line segment 
pq, called a feed-link, such that the largest network distance from p to any point in the resulting network is 
Q\ minimized. We solve the minimum eccentricity feed-link problem using eccentricity diagrams. In addition, we 

t-H provide a data structure for the query version, where the network G is fixed and a query consists of the point p. 

cn 1 Introduction 

J> We are given a network, i.e., an undirected graph with positive edge weights. We consider the continuum of 
k>( points along the edges of this network and measure distance between these points in terms of the shortest path 
$_J distance along the network. Within this metric space, we study farthest points. 

We introduce network farthest-point diagrams, which capture how the farthest points — and the distance to 
them — change as we traverse the network. We preprocess a network G such that, when given a query point 
q on G, we can quickly determine the farthest point (s) from q in G as well as the farthest distance from q in 
G. Furthermore, we introduce a data structure supporting queries for the parts of the network that are farther 
away from q than some threshold R > 0, where R is part of the query. 

This has applications to location analysis: Think of the network as roads in a city. When choosing the location 
for a new service facility, an urban engineer might want to know the parts of the network that are close-by 
(well-served) and the parts that are far away (ill-served). For instance, the farthest distance from a hospital 
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1 Introduction 





(a) The dilation of p and a point r on G is the ratio 
between the length of a shortest path (orange) from 
r to p via the extended network G + pq versus the 
Euclidean distance of p and r (blue, dotted). 



(b) The eccentricity of p in the extended network G + pq 
is the length of the feed-link (light red) plus the 
eccentricity of q in G, i.e., length of a shortest path 
(light green) from q to one of its farthest points q in G. 



Figure 1: An illustration of the minimum dilation feed-link problem (a), the minimum eccentricity feed- link 
problem (b), and their target functions. In both cases we have a network G and a point p in the 
plane that is connected to G via the feed-link pq (dashed) resulting in an extended network G + pq. 

influences the worst case response time of an emergency unit send from that hospital. Depending on the type 
of facility, we may be interested in locations with minimal farthest distance (network centers) or locations with 
maximal farthest distance (peripheral points) or points with farthest distance from some given location. 

We use the network farthest-point diagrams to solve the following network extension problem. Given a network 
G with geometric edge weights and a point p that is not on G, connect p to a point q on G with a straight line 
segment pq, called a feed-link, such that the largest distance from p to any point in the resulting network is 
minimized. In terms of our example, this corresponds to the task of connecting a hospital to a network of roads 
minimizing the worst-case emergency unit response time. The main difficulty of feed- link problems [4] stems 
from allowing feed- links to connect to any location in the network — not just to a few candidate locations — and 
taking every point on the network into account when measuring the utility of the feed-link. 



1.1 Related Work 

In their original work, Aronov et al. [4] introduced the feed- link problem with this example, i.e., connecting a new 
hospital to a network of roads. They consider a target function to measure the utility of a feed- link, other than 
its length. Aronov et al. [4] seek a feed-link that minimizes the worst-case detour one may take from any point 
on the network by traveling along the roads to the hospital as opposed to flying directly. The detour, or dilation, 
between two points p and q on a network is measured as the ratio of the distance of p and q via the network and 
the length of the straight line connecting p and q. We refer to this problem as the minimum dilation feed-link 
problem. Figure 1 illustrates dilation and the travel time to a farthest point (eccentricity) as target functions for 
the feed-link problem. See, for instance, Griine [15] for a comprehensive summary of dilation and its properties. 

We evaluate a feed-link with respect to all locations on a network. This means, if we use embedded graphs 
to model the network, all (uncountably many) points on this embedding — and not just (finitely many) vertices — 
count as possible farthest points. Aronov et al. [4] point out that the restriction to vertices yields a related 
feed-link problem where we seek the optimal feed-link to minimize the detour to a new train station in a railway 
system. In this scenario, we are only interested in the detour of the paths from other train stations to the new 
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one. If we model the stations as vertices, then the target function is the stretch factor [23] of the new station 
with respect to the extended railway system. 

Among other results, Aronov et al. [4] solve the minimum dilation feed-link problem for polygonal cycles 
in 0(A7(n)logre) time, where X-/(n) is the maximum length of a Davenport-Schinzel sequence [3] of order 
seven in n symbols. Agarwal, Sharir, and Shor [3] show that Xj(n) is almost linear in n, more precisely 
M( n ) < n • 2°( Q '™) loga(n))^ w j iere a ( n ) i s the inverse Ackermann function. Davenport-Schinzel sequences occur 
in the time bound, because Aronov et al. [4] rely on computing certain upper envelopes [2, 18]. Savic and 
Stojakovic [26] improved the result for polygonal cycles to 0(n) with their sliding lever algorithm. 

The subdivision of a network into parts with common farthest points is the farthest-point Voronoi diagram [25] 
on the metric space formed by the uncountably many points on the network and the network distance. Previously 
studied network Voronoi diagrams [11, 13, 16, 22, 24, 25, 29, 31] subdivide a network with respect to only finitely 
many reference points, e.g., depending on which reference point is closest or farthest. Okabe et al. [25] survey 
various notions of Voronoi diagrams, including some for networks. We refer the reader to Okabe et al. [24] for 
more types of network Voronoi diagrams and for further references. 

Creating maps of the farthest points in a network relates to problems from location analysis. For instance, 
in the continuous absolute 1-center problem [12] we seek a point in a network with minimal network distance 
to its farthest points. In terms of our prototype example, this is the problem of identifying the ideal position 
of a new hospital on a network. Kincaid [20] and Tansel [30] survey related notions and results. 

1.2 Preliminaries and Problem Definition 

A network is a simple, weighted, finite, connected, and undirected graph G = (V, E), where V is a finite set of 
points in R 2 , and E is a set of line segments whose endpoints are in V. Each edge e £ E has a positive weight 
w e > 0. We write uv for an edge with endpoints u and v. A point p on an edge uv G E subdivides uv into two 
sub-edges up and pv such that w up = \w uv and w pv = (1 — \)w uv for some A G [0, 1]. 

Consider the weighted shortest path distance dc ■ V x V — > [0, oo) between vertices of G with respect to 
the edge weights w e , e £ E. This can be extended to arbitrary points p and q on G by considering them to 
be vertices for the sake of evaluating dc(p,q)- More precisely, if we subdivide the edges containing p and q 
according to the above, then dc{p, q) is the weighted shortest path distance of p and q in the subdivided network. 
An example is shown in Figure 2. We refer to this distance as the network distance [4, 15] on G. 




u v 

(a) A network G. 




(b) Subdivision of G with respect to p and q. 



Figure 2: A network G (a) with edge weights as indicated. The subdivision of G according to the positions 
of p = \u + jv and q = hs + \t (b). We have dc{p, q) = 10 achieved on the highlighted path (blue). 
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Lemma 1. The network distance d G (-, •) is a metric on G, i.e., we have 

d G (p, q)>0 for all p,q £ G, 

d G (p, q) = if and only if p = q, 

d G (p,q) = d G (q,p) for all p,q G G, 

d G (p, q) < d G (p, r) + d G (r, q) for all p,q,r G G . 



(non-negativity) 

(identity of indiscernibles) 

(symmetry) 

(triangle inequality) 



Proof. Since all edge weights are positive, all paths on G have non-negative weight and only paths consisting 
of a single point have weight zero. Thus, non-negativity and the identity of indiscernibles hold, i.e., we have 
d G (p, q) > for all p,q G G where d G (p, q) = if and only if p = q . Since the edges are undirected, any path 
from p to q is also a path from q to p. Thus, we have symmetry, i.e., d G (p,q) = d G (q,p). Since d G (p,q) is the 
length of a shortest path from p to q, there cannot be a point r on G such that d G (p, q) > d G (p, r) + d G (r, q). 
Otherwise, the concatenation of a shortest path from p to r with a shortest path from r to q yields a path from 
p to q that is shorter than the shortest path from p to q. Thus, the triangle inequality holds. □ 

Let S C G be a subset of the points on G. The boundary of S, denoted by dS, consists of all points p G S 
such that for each value e > there exists one point on G within network distance e from p that is in S and 
there exists another point on G within network distance e from p that is not in S. For example, the boundary 
of the set of points on the highlighted path in Figure 2b consists of the points p, q, and x. 

The following definition generalizes graph eccentricity [17, pp. 35-36], which is usually introduced with respect 
to distances between vertices. Refer to Figure 3 for an illustration of this definition. 

Definition 2 (Eccentricity). Let G be a network and let p be a point on G. The largest network distance 
towards p is called the eccentricity of p with respect to G and it is denoted by ecc G {p), i.e., 

ecc G (p) := maxd G (p,q). 

A point q on G is called eccentric to p if it is farthest from p with respect to the network distance on G, i.e., 
if d G (p, q) = ecc G (p). We say that a point q is eccentric if there is a point p on G such that q is eccentric to p. 




(a) A network G. 




(b) A point p on G and a point p eccentric to p. 



2 u + 7}V with its eccentric 



Figure 3: A network G (a) with (geometric) edge weights as indicated. A point p — 

point p = f;S + ~t is shown (b). The eccentricity of p is ecc(p) = |\/5 + 4\/2 + 3, achieved on the 
highlighted path (orange). Neither p nor p are vertices of G in this example. 
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In the remainder of this article we will omit the subscript indicating the underlying network G in all of the 
above notation when it is clear from the context. Unless stated otherwise, we refer to the network distance 
whenever we describe distance between points on a network. 

Our terminology for ill-serviced parts of a network is defined as follows. Let p be a point on a network G, 
and let R > 0. We say a point q on G is R-far from p, if d(p, q) > R. Likewise, we call the set of points on 
G that are .R-far from p the i?-far sub-network of G from p. 

Alongside with our goal to characterize farthest-point information and changes therein, we aim to design 
efficient data structures to answer the following types of queries for any point p on a given network G. 

1. What is the eccentricity of pi 

2. Which points on G are farthest from p with respect to the network distance? 

3. What is the i?-far sub-network of G from p for some value R > 0? 

We represent the query point p as a pair made of an edge uv and a value A G [0, 1] with p G uv and p = (1— X)u+Xv. 
We refer to a query for the eccentricity as an eccentricity query, and we refer to a query for the set of farthest 
points as farthest-point-set query. 

With the notation established above, we can formally define the term feed-link and the minimum eccentricity 
feed-link problem as follows. An example is shown in Figure 4. 





(a) A network G with a point p ^ G. 



(b) An optimal feed-link. 



Figure 4: An instance of the minimum eccentricity feed-link problem (a). The anchor q of an optimal feed-link 
(dashed) (b) balances the distance to the eccentric points q and q' . 



Problem 3 (Feed-links and the Minimum Eccentricity Feed-Link Problem). Let p G R 2 be a point and 
G = (V, E) be a straight-line embedded network with geometric edge weights. 

(i) A straight-line segment pq with q G G is called a feed-link connecting p to G and q is called the anchor 
of this feed-link. The network that results from subdividing the edge containing q at q and adding the 
edge pq is denoted by G + pq. It is referred to as the extension of G by the feed-link pq. 

(ii) We call the task of finding a point q on G such that the feed- link pq with anchor q minimizes the eccentricity 
of p with respect to G + pq, i.e., the task of finding a point q on G that minimizes the expression 



ecc G+pq (p) = \pq\ + ecc G (q) = \pq\ + maxd G (g, r), 



(1) 



the minimum eccentricity feed-link problem. 
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1.3 Structure and Results 

The farthest-point Voronoi diagram of a set S of n points in K 2 subdivides the plane into regions with a common 
farthest point among those in S. Given some point p € R 2 , this diagram can be used to determine the farthest 
point from p in S . We wish to subdivide a network in a similar fashion in order to compute the farthest points 
from any point p on a network. However, there are differences between the situation in the plane and on a network. 
For instance, when given a network, the locations of the farthest points are unknown, as shown in Figure 5. 





(a) Stationary farthest points. 



(b) Moving farthest points. 



Figure 5: Two networks with different farthest point behaviour. In the network in (a), each point has one of 
the four squares as its farthest point. In the network in (b), each point p on the vertical middle edge 
has two farthest points p and p' . Moving p upwards causes these points to move downwards. 

Overcoming these issues and defining an analogue to the farthest-point Voronoi diagram in networks is one of 
the contributions of this work. We introduce the new notions of eccentricity diagrams in Section 2 and network 
farthest-point diagrams in Section 3. Network farthest-point diagrams encode the location of farthest points, 
whereas eccentricity diagrams capture the network distance to them. We show that, contrary to intuition, these 
diagrams may have non-linear size. In Section 4, we design and analyze a data structure for efficient eccentricity, 
i?-far, and farthest-point-set queries in networks. 

We solve the static and query version of the minimum eccentricity feed-link problem in Section 5 using 
eccentricity diagrams. We rephrase the query version as a point location problem in a certain Voronoi diagram 
whose sites are the sub-edges of the network with minimal eccentricity. Solving the static version of the feed-link 
problem takes 0(m 2 log n) time and 0(m 2 ) work space for a network with n vertices and m edges. Table 1 
summarizes the asymptotic bounds of all the other results. 



Query Time Pre-Processing Time 



Space 



Eccentricity Query 
i?-Far Query 
Farthest-Point-Set Query 



O(logn) 

0(fc + logn) 
0(fc' + logn) 



0(m 2 logn) 
0(m 2 logn) 
0(m 2 logn) 



0{m 2 ) 
0(m 2 logn) 
0(m 2 logn) 



Feed-Link Query 
Feed-Link Query 
Feed-Link Query 



0(£) 

o(io S e) 

OQogt) 



o(e 2+e ) 

expected 0{£ 2 log £) 



0(1) 
0{£ 2 ) 

o(e 2 ) 



Table 1: Our results for a network with n vertices, m edges. For R-iax queries, k is the number of edges 
containing i?-far points, and for farthest-point-set queries, k' is the number of edges containing farthest 
points. For the results regarding feed-link queries, £ S 0(m 2 ) is the number of sub-edges with locally 
minimal eccentricity, and we assume that these sub-edges are known a-priori. 
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In this section, we describe the distance to farthest points. We begin with computing the distance functions 
from points on an edge uv to their farthest points on an edge st. Combining these functions yields the distance 
from the points on uv to their farthest points in the entire network. This approach leads us to a representation 
of the distance to farthest points (eccentricity), which we call the eccentricity diagram of a network. From 
eccentricity diagrams, we derive a data structure for eccentricity queries, which we extend to a data structure 
for farthest-point-set queries in Section 3. 

2.1 The Shape of the Eccentricity Function 

We use a result by Frank [12] to compute the eccentricity of all points on a network G = (V,E). Frank [12] 
seeks a point on a network with smallest distance to its farthest points, i.e., a minimum of the eccentricity. 
He finds this minimum by determining a point with minimum eccentricity on each edge uv and then picking 
a point with the smallest eccentricity among these candidates. Frank [12] computes the eccentricity on an edge 
as follows. Let uv and st be edges of G. We define <f>^ v '■ [0, 1] — > [0, oo) as the mapping from A e [0, 1] to the 
largest network distance from the point p on uv with p = p(X) = (1 — X)u + \v to any point q on edge st, i.e., 

</>uvW := maxd(p(A),g) = maxd((l - X)u + Xv,q) . (2) 

qGst qGst 

The upper envelope of the functions 4>uvi over a ^ edges st e E, is the eccentricity of p(X) on uv, since 

ecc(p(X)) = max d(p(X) , q) = max max d(p(X),q) = max</>fL (A) . 

<jSG st&E qest st€E 

We begin the analysis of the functions (fi^ v an d their upper envelope with an auxiliary lemma. 

Lemma 4. Let G be a network, let ab be an edge of G, and let x be a point on G that is not in the interior 
of the edge ab. The network distance from x to the farthest point x from x among all points on ab is 

d(x x) - maxdix u) - d ( x ^ a ) + w ^ + d{x,b) 

a y X , X J — lllcLX. tiyx , y J — — . I O J 

ydab 2 

Proof. Let x be the farthest point from x on ab. Then the shortest path from x to x via a has the same length 
as the one via b, i.e., d(x, x) = d(x, a) + w a x = d(x, b) + Wbx- This yields 

_ 2d{x,x) d(x,a) + w ax + d(x,b) + Wxb 

d(x,x) = = , 

v ; 2 2 

and the result follows as w a b = w a s + Wxb- C 

The following lemma by Frank [12] describes the function <f>^ v for two distinct edges uv and st. Refer to 
Figure 6 for an illustration of the notation and the result. 

Lemma 5 (Frank [12]). Let uv and st be two distinct edges of a network G. Let u (respectively v) be the farthest 
point from u (respectively from v) on st. Likewise, let i (respectively s) be the farthest point from t (respectively 
from s) on uv. Without loss of generality, we have i G us (otherwise swap s and t). Then we have 



CXA) 
where p(X) = (1 — X)u + Xv and X £ [0, 1]. 



Xw uv + d(u,u) , if p(X) £ ut 

d(t,u) ,if P (X)eis, (4) 

(1- X)w uv + d(v,v) ,ifp(X)Gsv 
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d(v,t) 



i / / d(u,t) , 




d (t,u) = d(s,v) 



d(v,v) - 



\ _d(v,a) 



(a) Two edges in a network. 




d(u,u) - 



(b) The function p(X) h-> 4>^ v (X) — max gGst d(p(A), q), where 
p(A) = (1 - X)u + Xv and A e [0, 1]. 



Figure 6: An illustration of Lemma 5 and its notation for distinct edges uv and st. On the left (a), the edges 
are subdivided with respect to s, the farthest point from s on uv; to i, the farthest point from t on 
uv; to u, the farthest point from u on st; and to v, the farthest point from v on st. The positions 
of s, t, u, and v follow from Lemma 4. On the right (b), we see a plot of the network distance from 
points p on uv to their farthest point on st, according to Lemma 5. The slopes of the ascending and 
descending segments are equal up to their sign. 

The proof of Lemma 5 was omitted by Frank [12]. We add it for the sake of completeness. 
Proof. The three cases of Lemma 5 are illustrated in Figure 7. Case (3) is symmetric to Case (1). 



u = p 






(a) Case (1): p <G ut. 



(b) Case (2): p£ ts. 



(c) Case (3): p £ sv. 



Figure 7: An illustration of the three cases in Lemma 5 from left to right. The paths from p to p entering the 
edge st via s (orange) and via t (blue) are highlighted. 

Let p be the farthest point from p = p(X) on st. 

Case (1): Let p be located on ut as shown in Figure 7a. 

We show that u is the farthest point from p on st, i.e., we show p = u. Since p £ ut, there is a shortest 
path from p to t that includes the sub-edge up. Therefore, we have d(p,t) = w up + d{u,t). Likewise, 
we have d{p, s) = w up + d{u, s), since p is on us, as well. With Lemma 4 we obtain 



d{p,p) 



d(p, s) + d(p, t) + w st 
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w up + d(u, s) + w up + d(u, t) + Wgt 



w 



up 



+ 



d(u, s) + d(u, t) + w s t 



= w up + d(u,u). 

Case (2): Let p be located on is as shown in Figure 7b. 

We show d(p,p) = d(u, t) = d(s, v), which holds because p moves from u to v as p moves from t to s. 

Since p G iv and p G us, there is a shortest path from p to t that includes pv and there is a shortest 
path from p to s that includes pu. Therefore, we have d(p, t) = w pv + d(v, t) and d(p, s) = w up + d(u, s). 
Plugging this into (3) from Lemma 4 yields 



d(p,p) 



d(p, s) + d(p, t) + w st 

2 
w up + d(u, s) + w pv + d(v, t) + w st 



_ w uv + d(u, s) + ri(u, t) + w st 
2 

We know from Cases (1) and (3) that u is farthest from i on st, and that v is farthest from s on si. 
As the above applies to the cases p = i and p = s, we obtain d(p,p) = d(i, u) = d(s, v). 



Summarizing the three cases yields 

KiW =maxd(p(A),g) 



qest 



w up (\) + d(u, u) , if p(A) G ut 
d(t,u) , if p(A) 6 is , 

^w p (\)v + d(v,v) , if p(\)£sv 



where p(A) = (1 — A)n + Au, which implies the claim, since w up ^ = Xw uv and w p rx\ v = (1 — X)w v 



D 



In the next lemma we describe the distance from a point p on edge uv to its farthest point on uv itself — and 
thus the function </>™. An illustration of Lemma 6 and the notation used is shown in Figure 8. 

Lemma 6. Let uv be an edge of a network G. Let u (respectively v) be the farthest point from u (respectively 
from v) on uv. Further, let c be the midpoint of uv. Then we have u £ cv and v S uc with w U v = w vu , and 



C(A) 



' w uv + d(u, v) 

2 
(1 - \)w uv 

Xw uv 

w uv + d(u, v) 



if p{\) G uv 

if p(X) £ vc 

if p(X) G cu 

if p(\) G uv 



where p(X) = (1 — X)u + Xv and X G [0, 1]. 
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__ d(u, v) 



2 






U V 

(a) An edge in a network. 



(b) The function p(X) h* 0™(A) = max, eu „ d(p(X),q), where 
p(A) = (1 - A)u + Xv and A G [0, 1]. 



Figure 8: An illustration of Lemma 6 and its notation for an edge uv. On the left (a), the edge uv is subdivided 
with respect to u, the farthest point from u on uv; to c, the midpoint of uv; and to v, the farthest 
point from v on uv. On the right (b), we see a plot of the network distance from points on uv to 
their farthest point on uv itself, according to Lemma 6. The slopes of the descending and ascending 
segments are equal up to their sign, and the function 0™ is symmetric with respect to c. 



Proof. First, we show the claims about the positions of u and v. Using Lemma 4, we have 



»'(, 



,/ _ N Wuv + d(u, u) + d(u, v) w uv + d(u,v) w uv + d(v , v) + d(u, v ) ,, _. 
d(u,u) = = = = d(v,v) 



2 2 2 

which implies u G cv and v £ uc, since < d(u, v) < w uv . Furthermore, we have 



Wv 



U'„ 



w,, 



Wv 



WU1) Wq. 



Wf, 



Let p be the farthest point from p = p(X) on uv. 
Case (1): Let p be on uv with p/ti. 

We show that d(p,p) = d(u, u) = d(v, v) and that p moves from u to v as p moves from u to v. 

This case requires d(u,v) < w uv , as otherwise u = v. Let -k uv be a shortest path from u to v. The 
path ir uv and the edge uv form a simple cycle of length w uv + d(u, v). The farthest point p from p on 
this cycle has network distance d(p,p) = Wuv „ ^"'^^ = d(u,u) = d(v,v) to p. The points u, p, v, and 
u appear in this order along this cycle. Therefore, p appears between u and v, which shows p = p. 

Case (2): Let p be on vc. 

We show that v is farthest from p on uv, i.e., p = v. 

If we walk from u to p along uv, the distance to u increases from to w up , until we reach u. Hence, for 



all points q £ up we have d(p, q) < d(p, u) 



w 



up ■ 



If we walk from v to p along uv, the distance to v 



increases from to w vp , until we reach v. Hence, for all points q' £ pv we have d(p, q') < d(p, v) 
Since p S uc, we have w up < w pv and, thus, infer that p = v and d(p,p) 



w 



pu- 



ll) 



pv- 



The cases, p G cu and p G uv, are symmetric to Case (1) and (2), respectively, because w U v = Wu V . In summary, 

d(u,u) , if p(X) G uv 



^(A)=maxd(p(A),g) 



W up (\) 

d(y,v) 



, if p(A) G vc 
, if p(A) G cu 
, ifp(A) G nu 



which implies the claim, since w p ^x)v = (1 — A)w nt; , w up ^x) = \w uv , and d{u,u) = d{v,v) 



-d(u,v) 



D 
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2 Eccentricity Diagrams 

The eccentricity along an edge uv of a network with m edges is the upper envelope of m — 1 functions of the 
form described in Lemma 5 and Figure 6b and one function of the form described in Lemma 6 and Figure 8b. 
Thus, it is continuous and piece-wise linear. Next, we bound the number of linear pieces. To break ties, we 
number the functions (f) s ^ v and say that the function with higher number is higher wherever two functions coincide. 



Lemma 7. Let uv be an edge of a network G, and let k be the number of edges containing farthest points from 
some point on uv. The eccentricity on uv consists of 0(k) segments. 

Proof. For each non-constant segment only the part above the highest intersection with the other segments 
may appear on the upper envelope, due to the common slopes of the segments. Thus, each of the at most 
Ik non-constant segments contributes at most two bending points to the upper envelope. Since there is no 
intersection between two constant segments, this accounts for all bending points of the eccentricity function on 
uv, except for the first and the last one. Therefore, the eccentricity function on uv has at most 4k + 2 bending 
points and, thus, consists of at most 4k + 1 = 0(k) segments. □ 

2.2 The Eccentricity Diagram 

Due to the piece-wise linearity of the eccentricity, it suffices to state its value at the points corresponding to 
the endpoints of linear segments of the upper envelope of the functions (fi^v This describes the eccentricity 
on the entire network: For a point p with p = (1 — \)a + Aft on a sub-edge ab with linear eccentricity we have 
ecc(p) = (1 — A) ecc(a) + Aecc(6). This leads us to the following notion, which is illustrated in Figure 9. 

Definition 8 (Eccentricity Diagram). Let G be a network. We call the subdivision of G with linear eccentricity 
on every edge and with the minimum number of vertices the eccentricity diagram of G and denote it by £T>(G). 




A network. 



(b) Its eccentricity diagram. 



(c) Eccentricity values. 



Figure 9: From left (a) to right (c): A (geometric) network G, its eccentricity diagram £T>(G) with sub-edge 
weights, and its eccentricity diagram with the eccentricity at each vertex of £T>(G). 

The eccentricity diagram of a network is well-defined and unique, as it can be obtained by subdividing each 
edge uv at the finitely many endpoints of the line segments of the eccentricity function on uv. This yields a 
finite subdivision with the minimum number of additional vertices. An example is shown in Figure 10. 

Lemma 9. The eccentricity diagram of a network G = (V, E) with n vertices and m edges has size 0(m 2 ) and 
can be constructed in 0(m 2 logn) time. 

Proof. The eccentricity function along each edge uv consists of 0{k) = 0(m) segments, where k is the number 
of edges containing farthest points from uv. As we subdivide each of the m edges at the bending points of the 
eccentricity function, we obtain the bound on the size. 
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2 Eccentricity Diagrams 




(a) A network. 



e 2 



x 2 h 





1 2 




ft s V 

(b) The subdivision of uv. 



(c) The functions p(A) i->- (f>^ v (X) = max qest d(p(X),q), 
where p(A) = (1 — A)u + Av and A e [0, 1], together 
with their upper envelope. 



Figure 10: An example of how the edge uv of the network in (a) is subdivided (b) in the eccentricity diagram. 
In (c) we see a plot of the functions p(X) (->■ (p^ v (X) from points p on uv to the network distance 
to their farthest points on each edge st of the network. The upper envelope (thick, black) of these 
functions is the eccentricity on uv. We subdivide the edge uv at the points h, x, y, and a to achieve 
linear eccentricity on each sub-edge. 

We can construct the eccentricity diagram as follows. First, we compute the network distances between all 
pairs of vertices of G. We can use, for instance, Johnson's all-pairs shortest path algorithm [19], which has a 
running time of 0(n 2 log(ra) + nm) = 0(m 2 logn). With this information and Lemmas 4 to 6, we can determine 
the functions <p s ^ v for all edges uv,st £ E. Then we can compute the upper envelope of the functions (p^ v over 
all edges st £ E. For instance, Hershberger's Algorithm [18] can accomplish this task in O(mlogn) time per 
edge uv. The overall construction time is 0(m 2 logn). □ 

Next, we establish a lower bound on the size of eccentricity diagrams. The corresponding construction in 
the proof of Lemma 10 below shows that the bound stated in Lemma 9 is tight for planar networks. 



Lemma 10. For every k, I € 1ST with k > 2 and < / < k 2 , 
edges whose eccentricity diagram has size £l(kl). 



there is a network with 4k vertices and 4k — 1 + I 



Proof. Let k > 2 and let e be such that < e < ^— p Consider the network G/^o formed by the black edges 
of the network depicted in Figure 11. We obtain the network Gt.i by adding / edges of the form ViUj with 
i,j G {1, 2, . . . , k} and weight w ViU . = k — 1 to G^fl- All edge weights are positive, since k > 2. Thus, the 
network distance is a metric on Gt,i- 

Figure 12 illustrates the following arguments. First, consider a non-G/^o edge ViUj with i < j. We show that 
ViUj will be subdivided into at least 2k — 2 sub-edges in the eccentricity diagram of Gj-j, as shown in Figure 12a. 
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2 Eccentricity Diagrams 
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Figure 11: The black edges form the network G^q. Adding I orange edges ViUj with edge weight w Vi 
and i, j £ {1, 2, . . . , k} yields a network Gk,i- 
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(a) The farthest point x r from vertex x r on the edge ViUj 
with i < j for all r = 1, 2, . . . , k. 



(b) The farthest point p s t from p S UjUj on another 
non-G fe;0 edge u s v t . 




























































fc + 1 + e • i±i 
fc + 1 + e-i 

fc+i + e-^ 



Uj = X k Xk-l X 3 X 2 Xl 

(c) The distance from any point p on UjVi to the farthest among xi,Xz, . 



,a:fe. 



Figure 12: An illustration of the arguments in the proof of Lemma 10. Any non-G^ edge UjVi in network G^i 
is subdivided (a) into k edges in the eccentricity diagram of Gj-,i- We show that the distance from 
any point p on UjVi to the farthest point p s t on another non-G^o edge u s vt (b) is smaller than the 
distance to the farthest point among xi, X2, ■ ■ ■ , Xk (c). 

For each r = 1, 2, . . . ,k, let x r be the point on v^Uj that is farthest from x r . Then we have d(x r ,x r ) = k + l + e^-, 
for each r = 1, 2, . . . , k — 1 and Xk = Uj with d(xk, Xk) = d(xk,Uj) = k + 1 + ei. 

A plot of the mapping p i-» maxj =1 <i(p, x r ) from the points p on UjVi to their farthest point(s) among x\, X2, 
. . . , x r is shown in Figure 12c. We claim that this function coincides with the eccentricity function, i.e., that 
xi, X2, ■ ■ ■ , x r are the only eccentric points from the edge UjVf. We observe that the maximum distance to any 
of the vertices xi, X2, ■ ■ ■ , x r is always at least k + \ + e • ^ > k + e. 
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2 Eccentricity Diagrams 

On the other hand, all points on the network have distance at most k from p, except for points on edges 
incident to any of x±,X2, ■ ■ ■ , Xk- Let p s t be the farthest point from p on the edge u s Vt as shown in Figure 12b. 
Here we have d(p,p s t) < k because 

,/ _ N dip, u s ) + dip, v t ) + w UaVt 
dip, Pst) - 



2 




w pUj + e • \j - s\ + Wpvi +e-\i 


t\ -\- w UgVt 


2 




w Ua v t + WujVi + e • (\j -s\ + \i- 


~t\) 


2 


■■ 


k l + e- 2 




fe - 1 + e • (fc - 1) < k. 





Therefore, Figure 12c shows the eccentricity along ViUj. If i = j then this function consists of 2k — 2 segments and 
Ik — 1 otherwise. This is true for any of the I non-G^o edges. Hence, there are at least 4/c — 1 + Z • [2k — 2) £ Vi{lk) 
edges in the eccentricity diagram of Gkj- □ 

Corollary 11. For all n, m £ M wit/i m £ w(n) and n>8, there exists a network with n vertices and m edges 
whose eccentricity diagram has size Q(nm). 

Proof. Let n > 8 and m £ ui(n). Consider Gkfi from the proof of Lemma 10 with k = \j\- By adding 
I = m — n + 1 edges of the form ViUj, i,j = l,2,...,k, we obtain a network G^i with k = 0(n) and / = @(m). 
According to Lemma 10, the eccentricity diagram of Gk,i has size Q(kl) = Q(nm). Subdividing m mod 4 edges 
°f Gk i yields a network with the same property and exactly n vertices and m edges. □ 

Corollary 12. For all n £ IN with n > 8, there exists a planar network with n vertices whose eccentricity 
diagram has size Q(n 2 ). 

Proof. Let n > 8. Consider Gfc,o from the proof of Lemma 10 with k = [^J £ B(n). By adding all edges of the 
form ViUi for i = 1, 2, . . . , k, we obtain a planar network Gf-k- According to Lemma 10, the eccentricity diagram 
of this network has size Q(k 2 ) = Q(n 2 ). Subdividing n mod 4 edges of Gk : k yields a network with the same 
property and exactly n vertices. □ 

2.3 A Data Structure for Eccentricity Queries 

Assume we are given the eccentricity diagram £T>(G) of a network G as well as the eccentricity ecc(x) of each 
vertex x of £T>(G). Then we can answer queries for the eccentricity ecc(p) of any point p on G: It suffices to 
identify the edge ab of £T>(G) containing p, since 

ecc(p) = I 1 — ) ecc(a) -\ — ecc(6). 

V w ab J W a b 

Recall that eccentricity queries consist of a point p on G and the edge uv containing p. For each edge uv of 
G, we store the vertices x £ uv of £T>(G) (e.g., in an array or in a balanced binary search tree) sorted by the 
fraction ^^ at which x subdivides uv. Then we can find the sub-edge ab containing p with a binary search 
for ^^ in O(logn) time, as there are at most 0(m) = 0(n 2 ) vertices of £T>{G) on uv. 

Theorem 13. Given a network G with n vertices and m edges. There is a data structure that can be constructed 
in 0(m 2 logn) time and has size Oim 2 ) supporting queries for the eccentricity of any point p on G in O(logn) 
time, provided that the edge uv of G containing p is given. 
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3 Network Farthest-Point Diagrams 

3 Network Farthest-Point Diagrams 

Apart from the network distance towards farthest points, we are also interested in their location: We seek to 
query for the set of all farthest points from any point g on a network G. This suggests that we subdivide the 
network into parts with a common set of farthest points and then find the part containing g. In order to analyze 
this approach, we formally define this subdivision, which is the farthest-point Voronoi diagram whose metric 
space consists of all points on G and the network distance dc(-,-) and whose sites are all points on G. 

3.1 Farthest-Point Network Voronoi Link Diagrams 

Definition 14 (Farthest-Point Network Voronoi Link Diagram). Let G be a network. 

(i) Let g be a point on G. The set of points on G to whom g is a farthest point is denoted by Vf ar _ net (g), i.e., 

V far -„et(2) :={ P £G:Vg'£G: d(p,g') < d(p,g)}. 
We call Vf ar _ ne t(fl') the farthest-point network Voronoi link cell of g. 

(ii) We obtain the farthest-point network Voronoi link diagram of G by subdividing G at each boundary point of 
the non-empty farthest-point network Voronoi link cells, i.e., at all points in the set S := UoeG ^Har-net (#) • 

(iii) We say that a farthest-point network Voronoi link diagram is finite if and only if S is finite. 

Finite farthest-point network Voronoi link diagrams are, by definition, (finite) networks with finitely many 
vertices. Infinite farthest-point network Voronoi link diagrams, on the other hand, are infinite networks with 
infinitely many vertices and degenerate edges that are reduced to single points, as shown in Figure 13b. Tradition- 
ally, Voronoi diagrams are determined by a finite set of sites. The farthest-point Voronoi diagram [25, Section 3.3] 
subdivides the plane into regions with a common farthest point among finitely many points in the plane. The net- 
work Voronoi link diagram [25, Section 3.8] subdivides a network into parts with common closest vertices among 
finitely many vertices. The farthest-point network Voronoi link diagram differs from other Voronoi diagrams in at 
least two ways. First, we are unaware which points are the sites, i.e., which points on a network are farthest points. 
Second, there are infinitely many farthest points when the farthest-point network Voronoi link diagram is infinite, 
as depicted in Figure 13. We cannot immediately apply known methods to produce farthest-point network 
Voronoi link diagrams, and we cannot use these diagrams for farthest-point-set queries in the infinite case. 





(a) Finite diagram. 



(b) Infinite diagram. 



Figure 13: The farthest-point network Voronoi link diagrams for two networks. Parts that have a common 
farthest point (square) are indicated in colour. In the finite case (a), the network is subdivided into 
regions with a fixed farthest point. In the infinite case (b), we have a different behaviour on the 
vertical edges (black): When the point p is moved upwards, its two farthest points p and p' move 
downwards accordingly. No two points on this edge have a common farthest point. 
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3 Network Farthest-Point Diagrams 

Characterizing the finiteness of farthest-point network Voronoi diagrams reveals how we can avoid the 
difficulties of the infinite case. The following auxiliary lemma will help us with this characterization. 

Lemma 15. Given an edge uv of a network G and a point g on G. The set of points on uv that have g as a 
farthest point, i.e., the set Vf ar . ne t(g) n uv, is a (possibly empty) interval on uv. 

Proof. Assume, for the sake of a contradiction, that the statement is false. Then, there are points a and b on 
uv such that g is eccentric to a and b but not to some point c on ab. Let g' be a farthest point from c in G. 
Figure 14 shows a sketch of this (impossible) constellation. 



* u s a c b 

~ » — » C^— ^— ~9 



V 



Figure 14: The impossible constellation where the points on an edge uv that have g as farthest point consists 
of two disjoint intervals (thick, black). There would be two points a and b in these intervals such 
that the farthest point c of the midpoint c of ab is different from g. 

We first argue why g cannot be contained in ab. Lemma 6 implies that if g 6 uv and if g is farthest from 
both a and b, then we have either g = u or g = v. Recall that no two points on sub-edges with constant value of 
0™ have the same farthest points on uv, and that all points on the sub-edges of uv with increasing (respectively 
decreasing) value of §™ v have v (respectively u) as their farthest point on uv. We have g ^ ab in either case. 

Let g be the farthest point from g on ab. Without loss of generality, c is located on ag. As we walk from 
a to c along uv, the network distance to any point q in the network can change by at most w ac , i.e., for all q on 
G we have \d(c, q) — d(a, q)\ < w ac . The distance to g increases by exactly this amount as the network distance 
to g is increasing on the sub-edge ag, which includes the sub-edge ac, i.e, d(c,g) — d(a,g) = w ac . 

We assume that g is eccentric to a and b but not to c, i.e., d(c,g) < d(c,g'), d(a,g') < d(a,g) = ecc(a), and 
d(b,g') < d(b,g) = ecc(g) for some point g' on G. Therefore, the increase in the network distance to g' on ac 
must be strictly higher than the increase in network distance to g on ac, i.e., d(c, g') — d(a, g') > d(c, g) — d(a, g). 
This contradicts the assessment that the network distance to g already achieves the maximum possible increase 
of w ac along ac. Therefore, g' cannot exist and the claim follows. □ 

Theorem 16. Let G be a network. The farthest-point network Voronoi link diagram of G is finite if and only 
if there is no edge ab in the eccentricity diagram £T>(G) of G such that the eccentricity is constant on ab. 

Proof. Let the farthest-point network Voronoi link diagram of G be finite. First, we show that there are only 
finitely many points on G that are farthest from some point on G. We then infer that no sub-edge with constant 
eccentricity exists in G, and thus, no edge with constant eccentricity in £T>(G). 

Each of the finitely many boundary points of the farthest-point network Voronoi link cell has at most one 
farthest point per each edge of G. Since there is no change to the set of farthest points within the farthest-point 
network Voronoi link cells, there are only finitely many points on G that are farthest from some point on G. Let 
51) <?2j ■ ■ ■ ,9k be these points. As we walk along an edge uv of G, the network distance to gi strictly increases 
from u to gi, the farthest point from gi on uv, and then strictly decreases until v, as illustrated in Figure 15. 
Since the eccentricity on uv is the upper envelope of the network distances to the points gi,g2, ■ ■ • , 9ki the 
eccentricity cannot be locally constant on uv and there is no sub-edge of G with constant eccentricity. 
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3 Network Farthest-Point Diagrams 
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(a) Network 



" 94 92= 93 « = ffl 

(b) Distances and envelope. 



Figure 15: On the left (a), a sketch of a network with finitely many farthest points gi,g2,93, and 54. On the 
right (b), a plot of the network distances (coloured) to the points gi on edge uv. The upper envelope 
(black) of these finitely many non-constant functions cannot be locally constant on uv. 




u = p 



u a p {, v 

Figure 16: A sub-edge ab of an edge uv of G with strictly increasing eccentricity. Any shortest path (orange) 
from a point p £ ab to a farthest point p from p leaves uv via u. Thus, p is also a farthest point of u. 

Conversely, let there be no edge in the eccentricity diagram of G with constant eccentricity. Let ab be an edge 
of £T>(G) that is a sub-edge of edge uv of G. The eccentricity strictly increases or strictly decreases from a to b. 
Assume, without loss of generality, that the former is true as shown in Figure 16. Now, let p be a point on ab, and 
let p be a farthest point from p in G. Every shortest path from ptop leaves uv through u. Thus, any farthest point 
from p is also a farthest point from u. Therefore, there cannot be more non-empty farthest-point network Voronoi 
link cells than vertices of G. Due to Lemma 15, these have at most two boundary points in the interior of each edge. 
Hence, the number of boundary points is finite and so is the farthest-point network Voronoi link diagram. □ 

3.2 Network Farthest-Point Diagrams 

Theorem 16 shows that there are two ways in which farthest points change as we move along an edge of a network. 
First, the farthest points can remain the same. This happens if we move along a sub-edge with ascending or 
descending eccentricity. Second, the farthest points can all move staying at the same distance. This happens if 
we move along a sub-edge with constant eccentricity. In both cases, the edges containing farthest points change 
at most finitely often. Knowing the edges containing farthest points suffices to reconstruct the farthest points. 
We obtain a finite representation of the farthest-point network Voronoi link diagram by subdividing the edges 
of the eccentricity diagram depending on which edges of the network contain farthest points. 
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3 Network Farthest-Point Diagrams 

Definition 17 (Network Farthest-Point Diagram). Let G be a network. Consider the subdivisions of the 
eccentricity diagram £T>(G) of G such that the points in the interior of every edge uv of the subdivision have 
a common set of edges of G containing their farthest points, i.e., for all edges e of G and all p,q G uv\ {u, v} 



3p G e: d(p,p) = ecc(p) 



3q G e: d(q,q) = ecc(q). 



Among these subdivisions, we call the one with the least number of additional vertices the network farthest-point 
diagram of G and denote it by J-T>{G). 

Network farthest-point diagrams are well-defined and unique: We subdivide an edge of the eccentricity diagram 
wherever the set of edges containing farthest points changes. This occurs at most 0(m) times, i.e., when one 
of the functions <f)^ v joins with or departs from the eccentricity function, as shown in Figure 17. 




(a) A network. 





h x v f 

(b) The subdivision of uv. 



h x y f e a v 

(c) The functions <puv an d their upper envelope. 



Figure 17: An example for the network farthest-point diagram for the network in (a). The upper envelope of 
the functions (p^ v , which is shown in (c), reveals which edges contain farthest points. On the left 
(b) shows the network farthest-point diagram on uv. The farthest points from uv are located at the 
dot (ted segments) of matching colour: For instance, the sub-edge xy contains the farthest point for 
each point on the sub-edge xy, and vertex i is the farthest point for all points on the sub-edge ux. 

Storing the edges containing farthest points with every edge of the network farthest-point diagram yields a 
data structure for farthest-point-set queries of size 0(m 3 ), since there are 0(m) farthest points for each of the 
0(m 2 ) edges of FV(G). We answer a farthest-point-set query for a point p on edge uv as follows. First, we 
identify the edge oh in J r T>(G) containing p, using binary search. A list of the edges containing farthest points 
from p is stored with ab. For each edge st in this list, we compute the farthest point from p on si. This takes 
constant time per edge using Lemma 4. Thus, if p has k farthest points, we can report them in 0(k + logn) time. 

Theorem 18. Given a network G with n vertices and m edges. There is a data structure with size and 
construction time 0{m?) supporting farthest-point-set queries on G in 0(k + logn) time, where k is the number 
of reported farthest points. 
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4 A Data Structure for Eccentricity, R-Far, and Farthest-Point-Set Queries 

4 A Data Structure for Eccentricity, R-Far, and Farthest- Point-Set Queries 

When analyzing the location of a service facility in a network of roads, we may want to determine the part of 
the network that is ill-served, i.e., farther away from the facility than some critical threshold R. We introduce 
a data structure for i?-far queries, which consists of a point p on a network G, and a value R > 0. We seek the 
part of G with network distance at least R to the query point p. 

We develop a data structure for .R-far queries from a fixed edge uv and then build this data structure for 
every edge. To answer R-far queries from edge uv we perform two tasks: The first task is to compute the part 
of an edge st consisting of the i?-far points from query point p, i.e., the points q on st with d(p, q) > R. The 
second task is to identify those edges of the network that contain R-iai points without inspecting all edges. 

Consider a point p(X) on edge uv with p(X) = (1 — X)u + Xv for some A G [0, 1], and consider another edge 
st. Let p(X) be the farthest point from p(X) on edge st. The set of R-far points on st is the sub-edge of points 
q with w q p(^x) < 4>uvW ~ -^j since the distance to p decreases from p to v and from p to u. 

We rephrase the task to find all edges containing i?-far points as ray shooting problem. An edge st contains 
i?-far points from p(X) if and only if p(X) is R-far, i.e., if and only if d(p(X) , p(X)) = 4>u V {X) > R. Thus, we seek 
the functions (f>^ v , for all edges st, whose height at A is at least R, or, in other words, who are intersected by 
a vertical ray r that shoots upwards from the point (X,R), as shown in Figure 18. We solve this ray shooting 
problem separately for the line segments of the functions <j) s * v with a common slope (w uv , zero, or —w uv ). We 
use segment trees [5, 6] and exploit that line segments with a common slope have a vertical order. 




(a) A network. 





(b) The R-i&r portion of the network. 



(c) The functions 4>uvi f° r ai ^ edges st, and an 
upwards shooting ray with apex (A, R). 



Figure 18: An example for the correspondence between an R-iar query in the network and a vertical ray stabbing 
query. On the upper left (a), we see a network G. On the lower left (b), we see the answer (red) 
to an i?-far query at p(X) = (1 — A)u + Xv on edge uv of G with A = j| and R = 13.5. On the right 
(c), we see the corresponding ray stabbing query. The functions <p^ v that are stabbed by the ray 
(red) contain R-far points from the query point p(X). 
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5 The Minimum Eccentricity Feed-Link Problem 

Lemma 19 ([5, 6]). Let L be a set of N vertically ordered line segments in K 2 . There is a data structure that 
reports all k line segments in L that are intersected by a vertical ray shooting upwards from a query point (x, y) 
in 0{k + logiV) time. This data structure has construction time and size 0(N log N). 

Our data structure for R-iar queries on an edge uv consists of three segment trees T asc , T con , and Td es containing 
the ascending, constant, and descending line segments of the functions (j)^ v for all edges st in decreasing vertical 
order. These segment trees allow us to compute the k edges of the network containing i?-far points from a query 
point p(X) = (1 — X)u + Xv in 0{k+\ogn) time with a ray shooting query for the ray shooting upwards from (A, R). 

Lemma 20. Let uv be an edge in a network with n vertices and m edges. There is a data structure that supports 
R-far queries from points on uv in 0(k + logn) time, where k is the number of edges containing R-far points 
from the query point. This data structure has size and construction time O(mlogn). 

We can use this data structure for eccentricity queries and farthest-point set queries on edge uv, as well. For 
eccentricity queries, we keep track of the maximum heights at A of the line segments stored at the heads of the 
lists encountered as we follow the search path for A. The maximum height is the eccentricity of p(X), since it 
is the greatest value among </>„„(A) for all edges st of G. This query takes O(logn) time, since we only have 
to inspect the heads of the O(logn) lists along the search path. For farthest-point set queries, we first determine 
the eccentricity of p(X) and then perform a R-far query with R = ecc(p(X)). We obtain our final data structure 
for all three types of queries by building the data structure from Lemma 20 for each edge of the network. 

Theorem 21. Given a network G with n vertices and m edges. There is a data structure with size and 
construction time 0(m 2 log n) supporting eccentricity queries, R-far queries, and farthest-point-set queries from 
any query point p on G. Let k denote the number of edges containing R-far points from p, and let kl denote 
the number of farthest points from p in G. Using the data structure, an eccentricity query takes O(logn) time, 
an R-far query takes 0(k + logn) time, and a farthest-point- set query takes 0(k' + logn) time. 

5 The Minimum Eccentricity Feed-Link Problem 

Ln this section, we solve the feed-link problem. 

Problem 3 (Feed-links and the Minimum Eccentricity Feed-Link Problem). Let p G R 2 be a point and 
G = (V, E) be a straight-line embedded network with geometric edge weights. 

(i) A straight-line segment pq with q £ G is called a feed-link connecting p to G and q is called the anchor 
of this feed-link. The network that results from subdividing the edge containing q at q and adding the 
edge pq is denoted by G + pq. It is referred to as the extension of G by the feed-link pq. 

(ii) We call the task of finding a point q on G such that the feed- link pq with anchor q minimizes the eccentricity 
of p with respect to G + pq, i.e., the task of finding a point q on G that minimizes the expression 

ecc G+pq (p) = \pq\ + ecc G (q) = \pq\ + max d G (q, r), (1) 

r&G 

the minimum eccentricity feed-link problem. 

We treat two versions of the minimum eccentricity feed-link problem. In the static version, we have a fixed 
network G and a fixed point p that we wish to connect to G. Here, we seek the optimal feed-link for G and 
p. In the query version of the problem, we have a fixed network G and a query consists of a point p that we 
wish to connect to G. Here, we seek a data structure that can answer queries of this type efficiently. 
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5 The Minimum Eccentricity Feed-Link Problem 



5.1 The Static Version 

The dependence upon the eccentricity ecc G (q) of the anchor point q in (1) shows the connection between farthest- 
point information on networks and the minimal eccentricity feed-link problem. This dependence determines 
necessary conditions on the optimal feed-link. 

Lemma 22. Let G be a geometric network. Let the point q on G be the anchor of an optimal feed-link for the 
point p E R 2 . Then the eccentricity on G has a local minimum at q. Furthermore, if q is located on an edge 
uv of £T>{G) with constant eccentricity, then q is the closest point to p on uv with respect to Euclidean distance. 

Proof. Let uv be the edge of £T>(G) containing the optimal anchor q of a feed-link from p to G. 
Case (1): Let the eccentricity be increasing on uv with ecc G (u) < ecC(j(v). Then we have 

ecc G+pg (p) = \pq\ +ecc G (<?) 

> \pu\ - \uq\ + ecc G (q) 

= \pu\ — \uq\ + ecc G (u) + \uq\ 

= \pu\ + ecc G (u) 

= ecc G + P u(p) ■ 



as the eccentricity is increasing on uv, 



Therefore, the optimal anchor among all points on uv is u, i.e., q = u. 

Case (2): Let the eccentricity be constant on uv, and let q' be the closest point from p on uv. 



ecc G+pq (p) = \pq\ +ecc G (q) 
= \pq\ +ecc G (q) 
> \pq'\ + ecc G (q) 
= ecc G+pg ' (p) . 



as the eccentricity is constant on uv, 
by choice of q' , 



Therefore, the optimal anchor among all points on uv is q' , i.e., q 



q'- 



□ 



Using Lemma 22, we can solve the minimum eccentricity feed-link problem as follows. First, we compute 
the eccentricity diagram £T>{G) of the network G. Second, we read the sub-edges of G with locally minimal 
eccentricity from £D(G). These sub-edges are the edges of £T>(G) with constant eccentricity and the vertices 
of £V(G) — which we treat as sub-edges reduced to single points — whose neighbours in £V{G) have greater 
eccentricity. Let S be the set of sub-edges with locally minimal eccentricity. Third, we determine a candidate q 
for the optimal anchor on each sub-edge with locally minimal eccentricity. Among these candidates we select the 
one with the lowest value of \pq\ + ecc G (q). An example is show in Figure 19. The first step takes 0(m 2 logn) 
time as discussed in Section 2. The second step can be done alongside with the computation of the eccentricity 
diagram. The third step takes 0(|S|) time, where \S\ S 0(m 2 ). 

Theorem 23. Given a geometric network G with n vertices and m edges, and a point p £ M 2 . Assume we 
are given the I sub-edges of G with locally minimal eccentricity. Then we can solve the minimum eccentricity 
feed-link problem with respect to G and p in 0(£) time. 



Corollary 24. Given a geometric network G with n vertices and m edges, and a point p G R 2 . 
the minimum eccentricity feed-link problem with respect to G and p in 0(m 2 logn) time. 



We can solve 
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(a) A network G with a site p. 



(b) Canditate feed-links. 



(c) The optimal feed-link. 



Figure 19: From left (a) to right (c): An instance of the minimum eccentricity feed-link problem. The candidates 
for optimal feed-links (dashed) and local minima (thick, black) of the eccentricity function. The 
subdivision of the plane into regions with a common sub-edge containing an optimal feed-link. 

5.2 The Query Version 

Now we address the query version, where the network G is fixed and a query consists of the point p. Throughout 
the following let S be the set of sub-edges of G with locally minimal eccentricity, and let £ = \S\. 

Using the solution for the static problem, we can create a data structure with size 0(1) and construction time 
0(m 2 logra) that answers queries for the optimal feed-link in 0{€) time. This data structure consists of the set 
S, and we obtain it by computing the eccentricity diagram and recording the local minima of the eccentricity. 
We improve the query time — at the expense of space consumption and construction time — by rephrasing the 
minimum eccentricity feed-link problem as a point location problem in a special type of Voronoi diagram. 

We denote the Euclidean distance between a point p G IR 2 and a segment s G S by d 2 (p, s), i.e., d 2 (p, s) = 
min ggs |pg|. By definition, the eccentricity with respect to the network G is constant on all segments s G S. We 
write ecc(s) to denote the eccentricity of the points on s. With this notation, the optimal feed-link is the closest 
sub-edge s £ S with respect to the additively weighted Euclidean distance d 2 (p, s) + ecc(s). 

Conversely, consider the Voronoi diagram of the line segments in S with respect to the additively weighted 
Euclidean distance where the weight of a segment s G S is its eccentricity ecc(s). This diagram splits the plane 
into regions whose points have a common closest segment in S with respect to the additively weighted Euclidean 
distance. In other words, the points in each region have their anchor of an optimal feed-link on a common 
sub-edge in S. Figure 19c shows an example of this kind of Voronoi diagram. 

Definition 25. Let S be a set of line segments in the plane with weights w s G H for each s G S. 

(i) The additively weighted distance of a point p £ M 2 and a line segment s G S is the Euclidean distance 
d 2 (p, s) of p and s plus the weight w s of the line segment s. 

(ii) We call the set of points p G R 2 to whom a line segment s G S has the lowest additively weighted distance 
among all line segments in S the additive weight Voronoi cell of s, and denote it by V+(s), i.e., 

V+(s) :={pGR 2 : Vs' G 5: d 2 (p,s) + w s < d 2 (p,s') + w s ,}. 
(iii) We call the subdivision of the plane into the set [J s£S dV+(s) and the connected regions of 



Vses / 



the additively weighted Voronoi diagram of the line segments in S with respect to the weights w s , s G S. 
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(a) Two line segments and their weights. 



(b) A part of the subdivision of the plane into regions 
(green) with smaller additively weighted Euclidean 
distance to ab than to cd and their complement (red) . 



Figure 20: The bisector (b) of two line segments with additive weights (a). The visible part of this bisector 
consists of line segments, parabolic arcs, and hyperbolic arcs. The region of points closer to the 
green segment are non-convex and disconnected; this misbehaving bisector splits the plane in three. 

Voronoi diagrams of points with additive weights [25, Section 3.1.2] and Voronoi diagrams of line segments 
[25, Section 3.5] have received considerable attention in the literature and are thus well studied concepts. To 
the best of the authors' knowledge, the following quotation is the only direct mentioning: 

"In general, the Voronoi diagram of segments where each segment carries an additive weight is not 
a well-behaved Voronoi diagram: Voronoi regions can be disconnected, and the diagram can have 
quadratic complexity." — Cheong et al. [8] 

The comprehensive study of additively weighted Voronoi diagrams of line segments is beyond the scope of 
this work. Nonetheless, we summarize a few observations about this type of Voronoi diagram. Figure 20 shows 
the (ill-behaved) additively weighted Voronoi diagram of two line segments. 

Theorem 26. The additively weighted Voronoi diagram of £ line segments is an planar subdivision of size Q(£ 2 ) 
whose edges are parts of lines (lines, rays, line segments), parts of parabolas (parabolas, parabolic rays, parabolic 
arcs), and parts of hyperbolas (hyperbolas, hyperbolic rays, hyperbolic arcs). 

Using common techniques from planar point location [9, 28], we determine the region that contains a query 
point p in 0(log£) time. The storage requirement and construction time of this data structure are linear in 
the size of the planar subdivision, provided that the subdivision is monotone. We can make the subdivision 
monotone in 0(£ 2 logi) time using plane sweep. 

We briefly discuss the construction of additively weighted Voronoi diagrams of line segments. The algorithms 
for abstract Voronoi diagrams [21] only work for (generalizations of) Voronoi diagrams whose Voronoi cells are 
connected, a property that is violated in our case. Despite the lack of existing theory, we can compute the desired 
diagram using the relationship between Voronoi diagrams and lower envelopes in three dimensions via lifting 
maps. Setter, Sharir, and Halperin [27] provide a comprehensive review of this observation by Edelsbrunner 
and Seidel [10] and its practical implications. In a nutshell the idea is as follows: consider the lower envelope of 
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the graphs of the distance functions of all sites of the Voronoi diagram. The projection of this envelope onto the 
plane yields the desired Voronoi diagram. Agarwal, Schwarzkopf, and Sharir [1] provide a divide-and-conquer 
algorithm that computes the lower envelope of the (weighted) distance functions of £ sites in 0(£ 2+e ) time. The 
randomized version of this algorithm, which was proposed by Setter, Sharir, and Halperin [27], accomplishes 
the same task in 0{£ 2 log^) expected time. 

Theorem 27. Let G be a geometric network with n vertices and m edges. Furthermore, let t be the number 
of sub-edges of G with locally minimal eccentricity. There is a data structure that can perform queries for a 
minimum eccentricity feed-link for any point p G IR 2 in 0(\og£) time. This data structure has a space requirement 
of 0(£ 2 ). It can be constructed in 0(£ 2+e ) time or, alternatively, in 0(£ 2 log£) expected time, both provided that 
the eccentricity diagram of G is known a-priori. 

6 Conclusions and Future Work 

We introduced new notions to capture farthest-point information in networks as well as data structures to store 
and access this information efficiently. We seek to improve the bounds on the construction time and space 
requirements in future work. For instance, our approach ignores any structure that the network might have 
and requires all pairs shortest path distances. 

We presented the feed-link problem that kindled this research alongside with a first solution for its static and 
query version. The feed-link problem can be extended in many ways. For instance, we could connect several 
sites simultaneously to a network minimizing the largest distance to the nearest site. Further, we could require 
the extension of a planar network to be planar as well or add other restrictions such as obstacles. Aronov et al. 
[4] discuss the latter two for the minimum dilation feed-link problem. 

Finally, the additively weighted Voronoi diagram of line segments demands more investigation, because of 
its relation to the query version of the feed-link problem. 
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